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^\ , Large-angle fluctuations in the cosmic X-ray background are investigated by a new 

0^ ' formalism with a simple model of the X-ray sources. Our method is formulated from 

the Boltzmann equation and a simple extension of the work by Lahav et al. to be 
applicable to a hyperbolic (open) universe. The low multipole fluctuations due to the 
' source clustering are analyzed in various cosmological models in both numerical and 

analytic way. The fluctuations strongly depend on the X-ray sources evolution model, as 
pointed out previously. It turns out that the nearby (z ^ 0.1) sources are the dominant 
CO ' contributors to the large-angle fluctuations. If these nearby sources are removed in an 

observed X-ray map, the dipole (low multipole) moment of the fluctuation drastically 
decreases. In this case the Compton-Getting effect of an observer's motion can be a 
dominant contribution to the dipole fluctuation. This feature of fluctuation, relating to 
the matter power spectrum, is discussed. 
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oo '. I- INTRODUCTION 

. One of the most important problems in cosmology is to explain the formation of the cosmic struc- 

CIh' ture. The discovery of anisotropics in the cosmic microwave background (CMB) supports the hypothesis 

that the structure of the Universe formed via gravitational instability from small density fluctuations 
■ (e.g., White, Scott, & Silk 1994). It is also argued that the cosmological models with cold dark matter 

^ I (CDM) are favorable models with which to explain the large-scale structure of the galaxy distribution 

d • and the formation of the cosmic objects at the high-z Universe (e.g., Dodelson, Gates, & Turner 1996; 

Peacock & Dodds 1994). 

The origin of the cosmic X-ray background (CXB) radiation is also explored (e.g., Fabian & Barcons 
1992;Hayashida 1989). The deep surveys of the 0.5-2 keV X-ray range were carried out with ROSAT, and 
most of the background sources were resolved into AGN and other X-ray luminous galaxies (e.g., Hasinger 
1997). From the isotropy of the observation, it was concluded that a good deal of the unresolved fraction 
of the CXB has come from redshifts z > 1 (e.g., Carrera, Fabian & Barcons 1997). The conclusion may 
be supported in the higher X-ray energies band with ASCA observations (Georgantopoulos et al. 1997). 
The finite solution of the source problem will be obtained in the next generation of X-ray instruments 
(e.g., Charles & Seward 1995). 

The CXB fluctuations are an important probe of the large-scale structure of the universe ( Barcons 
& Fabian 1988; Carrera et al. 1997; Barcons, Fabian & Carrera 1997; Lahav, Piran & Treyer 1997). 
For example, the relation between the angular correlation function of the intensity fluctuations and the 
two-point spatial correlation function of galaxies ^(r) was studied. (De Zotti et al. 1990; Martin-Mirones 
et al. 1991; Persic et al. 1989; Shafcr 1983). The cosmological models with a cosmological constant can 
be tested by the cross-correlation between the X-ray and the microwave background(Boughn, Crittenden 
k Turok 1997). 

Recently, Lahav et al.(1997) have investigated the large-angle fluctuations in the CXB. They have 
calculated the expected values of the fluctuations in a statistical way. Treyer et al.(1998) have compared 
the predicted CXB fluctuations with the HEA01-A2 measurements. In the prediction of CXB fluctuations 
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there are many uncertainties. The evolution of the X-ray sources are not completely understood. A simple 
(power-law) source evolution model is assumed in their paper (Lahav et al. 1997; Treyer et al.l998). 
Furthermore only a simple flat cosmological model is assumed. The difference in the cosmological model 
might yield the large difference in amplitude of the fluctuation. For example, if the X-ray sources at the 
high-z universe are the dominant contributors to the angular fluctuations, the cosmological parameters, 
e.g., the curvature of the universe and the cosmological constant A, are significant factors. It is therefore 
worth examining how the fluctuations depend on the cosmological parameters. We also develop a useful 
formalism which is applicable to a hyperbolic (open) universe in order to extend the work by Lahav 
et al. to various cosmological models, which is described in Section 2. Our formalism is based on the 
Boltzmann approach, which is familiar in analysis of the cosmic microwave anisotropies (e.g., Hu & 
Sugiyama 1995a;1995b). A simple model for the source distribution is introduced. In Section 3 we solve 
the Boltzmann equation and obtain the expression for the root mean square of multipole moments of 
the fluctuations. In Section 4 the fluctuations are analyzed in various cosmological models. Section 5 is 
devoted to summary and discussions. 

Throughout this paper we will use the units c ^ Ti = ks = I', however, we occasionally use the Planck 
constant hp{= 2'Kh) to make clear the meaning of equations. 



II. FORMALISM 



The perturbation of the Fricdmann-Robertson- Walker (FRW) space-time in the Newtonian gauge is 
written as 

ds^ = -(1 + 2^)dt^ + a^(l + 2<^)-iydx'dx' , (2.1) 

where ^ is the perturbed gravitational potential, $ is the curvature perturbation, 7.^ is the three-metric 
on a space of constant negative curvature, and a is the scale factor normalized to be unity at present. 
The scale factor a is given by the Friedmann equation, 

=H^a(^^o + ^K + a^nf}j , (2.2) 

where is the Hubble parameter with i/o = 100 h km/s/Mpc, fio is the density parameter, ilA(= 
A/3Hq) is defined by the cosmological constant A, and flxi— ^ — ^q — ^a) describ es th e spatial curvature 



of the universe. For the spatially flat models, flK is equal to zero. In equation (2_^) the dot denotes rj 
differentiation, where 77 is the conformal time defined by adrj — dt. 

The propagation of photon is affected by the expansion and the metric perturbations. The Boltzmann 
equation in the perturbed FRW universe is derived in Appendix A. The explicit form is expressed as 

dT]~^^ ^ ^ dx"^^ dq^a^ drj ^ dx''^ ) ^ d-f^ dr] iVi tIjI ] ^ ( ) 

where q is the photon energy in the locally orthonormal frame, and 7' is the spatial vector normalized 
as 7^7*7"' = 1 in this frame. 



We now discuss the source term S in the right-hand side in equation (2.3). The X-ray sources may be 
the discrete point sources. As we are interested in the much larger scale, we can average the distribution 
of the X-ray sources over a certain region and regard the source function as a smoothed field with spatially 
varying emissivity. We write the comoving volume emissivity per frequency at the time 77 as j^iVi ^) 
separate it into two parts as 

j^(77,x) =>(??) +Sji,{T],x.) , (2.4) 

where the bar represents a spatially averaged quantity, and Sj^, is the spatially inhomogeneous part of 
ji,. In this paper, we consider the cosmological models dominated by CDM. The emissivity may be 
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changed due to the presence of the inhomogeneity of the CDM. We expect that the value ji, should 
increase as the CDM density fluctuation 5c (= [pc ~ Pc)/Pc) becomes large. Therefore, the emissivity in 
the inhomogeneous universe is expressed as 

Mri,^)=Mrj){l + bx5c) , (2.5) 

where bx is the bias factor to relate the large-scale distribution of the CDM with that of the X-ray 
sources. |^ In this paper we assume 

J.iv) = Jo{'^)E{z) , (2.6) 

with E{z) = {l + z)P for redshift Zmin < z < z^ax, and E{z) — for other period. Note that p is constant 
and p = is the case of no evolution of the X-ray sources. We also assume the power-law form of the 
energy spectrum, jo{h') oc with a = 0.4. Here we set a to be the observed CXB spectral index in 
(3-20) keV range (Comastri et al. 1995), however, a reasonable change in a will not affect our conclusions 
at a significant level. 

The photon number emitted during the proper time t t + dt, in the frequency range v ^ + dv, in 
the solid angle dVt, from the volume x ~ x -f dx, is written as [J 



f 1 + (3 + a)-i^V'x) dv'^d^^dt , (2.7) 
/ipi^ \ / 47r 



where we assumed that this volume element has the peculiar velocity Vx^ and note that the Doppler 
effect is taken into account. 

From the definition of the source term S in the Boltzmann equation, the photon number emitted during 
t t + dt^ in the momentum range q ~ q-|-dq, from the volume x ~ x-|-c?x, is Sgd^qa^d^xdt/hp^ , where 
g{= 2) is the photon statistical weight. Using the relation q{= |q|) = hpiy, the photon number emitted 
during t t + dt, in the frequency range h' ^ ly + dv, from the volume x ^ x -I- dx, in the solid angle dU,, 



is written as gSv'^di'dila^d'^xdt. Equating this expression and equation (2/7), we find 

1 1 J. 



■(l + (3 + a)7.Vl), (2.. 



where we used the relation hp = 2tt in our unit. The peculiar velocity of the X-ray sources Vx does not 
necessarily agree with that of the CDM, Vc- We introduce another bias factor by for velocity, where Vc 
is the CDM velocity field. 

We eventually get the source term as 

S=^^^U + bxSciv, x) + {3 + a)bvj,Vc\fi, x)) . (2.9) 



* It may be unrealistic to assume the epoch-independent biasing. The epoch-dependent biasing is assumed in 
the recent paper by Treyer et al.(1998) by setting bx{z) = bx{0) + z[bx{0) — 1] . As we describe in section 4 and 
Appendix D, however, the nearby sources {z ^ 0.1) are the dominant contributors to the large-angle fluctuations 
(Z ^ 10). In this case the amplitude of large-angle fluctuations will roughly scale in proportion to bx{0) unless 
bx{0) becomes very large. 

^ The fluctuation of the metric might have to be taken into account when formulating of the source term. As 
we will see in the next section (also Lahav et al. 1997), however, the contribution from the metric fluctuations to 
the large-angle fluctuations in the CXB is small. Therefore we have omitted the metric fluctuations in the source 
term. 
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III. SOLUTION OF THE BOLTZMANN EQUATION 



In order to solve the above Boltzmann equation (3), we use the hnear perturbative expansion method 
with respect to the inhomogenity (such as metric perturbations) . If the universe is completely homoge- 
neous, the distribution function /'^"^ and the source term can be described by rj and q. With the 
inhomogenity, they are modified as 

f^f^'\v,q) + f'^'\v,^\q,ii , (3.1) 

5 = 5(°)(7y,g) + 5W(77,x\g,7'0 , (3-2) 

with, 

5(")(^,,) = 4^1:^, (3.3) 
5(1) (77, x\q, y ) = 5(") (ry, q) (bx5,{rj, x') + (3 + a)bv^^V,\f], x')) . (3.4) 
The Boltzmann equations of the zeroth and first order are given by 



drj 
and 

drj dx^ ^ dq a dq \ d?] dx^ ^ ) dj'^ drj 

Note that when neglecting the metric perturbations, equations (3.5) and ( |3.6| ) are equivalent to the 
original Boltzmann equation, and give the exact solution. 

Introducing the comoving momentum by Q = aq, the above equations reduce to 

^-a5(0)(r;,Q/a) , (3.7) 

and 

+ i^r-r - Q^^\ IT + TTii + -^r-r-r- = aS^ 'iv,x\Q/a,j') , (3.8) 



dri dx^ dQ \ drj dx'- J dY drj 

where /(°) and /(^^ must be regarded as the functions f''^\'r],Q) and /(i)(?7, x*, Qj 7*)j respectively. 
Let us first consider the zeroth order equation (|3.7D, which is formally integrated as 



Jo 

= -Q-' r d7^'a'jQ/^^,,{rj') , (3.9) 
9 Jo 

where we used S'--^\Q/a) ~ {t^ / gQ'^)3Q/h-pa{fi)j ^iid defined a' = a{ri'). Since we have assumed the 

power-law form, ji/(ry) = + zY , for Zmin < z < Zmax, (see equation (p^)), the specific intensity, 
/(,") ^ ghpv^f^^'i (e.g., Shu 1991) is written as 

= ^ r '^<'+^y-'- , (3.10) 



4 



where we used the relation be twee n the conformal time and the redshift (2.2). 

We next consider equation (3.8). It is convenient to introduce a new variable 0(77, x*, Q, 7') as 



(3.11) 



where we used the general relation /j, — ghpv^ f for the second equality. By using 8(77, a;*, Q, 7*), equation 
(|3.8[) can be written as 



drj 



de 



+ e(^)0 + TrT7' + (a + 3) irT7^ + ^ + 



dx 



9* 



de dY 



dx^ drj J dY dr] 
e(r7)(fex<5c(x) + 6y(3 + a)7,Fc*(x)) , 



(3.12) 



where we used the relation QOf'^^^dQ = —{a + 3)/*^°^ and equation (3/7), and ^ is defined by ^(r/) 

(a/(")/aQ)/(a/(0V5Q) = /("V/*°^ = eW/^W with 



e(ry) = / dij' ajQ/h^a'iv') / d-q' E{-q')a 



(3.13) 



Note that the variable Q is separated in equation (3.12) and that can be regard as the function 
0(77, X, 7'). This is the result of the simple assumption of the power-law energy spectrum ( ^.6| ). 
Equation ( |3.12 ) can be written as 



d_ 

drj 



;(e(ry, x(?7), 7(7/)) + (a + 3)^(77, x(77))^ 

= 2{a + 2,)e^/^^[e^/H{'q, x(7;))) +e(bx6,{i^, x) + 6y (3 + a)-i,V^) , 



(3.14) 



where we used $ = — ^' from the linear perturbation theory. This can be easily integrated as 
e(77o,x(?7o),7) + (a + 3)*(?7o,x(7/o)) 



drj 



(a + 3)£i/2^ (ei/2*(77, x)) +i{bx6civ, x) + 6y (3 + 0)7. K 



(3.15) 



where we used £(0) = 0. 

We rewrite this solution by using the multipole expansion method. The formalism for the CMB 
anisotropies is useful (e.g., Hu & Sugiyama 1995b). Since 9(77, x, 7*) follows the linear equation, we can 
expand it as a sum of modes labeled by (k, I), 



e(7y,x,f)=^^eKr/,fc)GKx,y) , 



(3.16) 



k ; 



where Gi is a certain mode function. Th e explicit form is written in Appendix B. The solution for 
ei{r]o,k) can be found from equation ( 3.15 ): 



Qoivo^k) = -(a + 3)*(7yo,A:) + 



1 



d77 



(« + 3)£V2^(ei/2vi.(^,fc))xO(A7?) 



sivo) Jo 



(3.17) 



and 



5 



1 



drj 



21 + 1 e(r/o. 



I 



with Arj 



21 + 1 
K{rjQ — rj) and 



i + i 



21 



l{l + 2)K-^,^, 



fc2 



a > 1) , 



2 sinh X 



1/ 



(3.18) 



(3.19) 



where Py{x) is the Legendrc function, and u is defined as = k"^ /{—K) — 1 with the eigen- value 
fc^ of scalar ha rmon ics an d the spatial curvature parameter K = —Hq{1 — ^Iq — JIa). We can check 
that equations ( 3.17| ) and ( 3.18| ) satisfy the perturbation equations for the multipolc moments derived in 
Appendix B. 

In the limit of flat u nivers e, i.e. K ^ 0, the radial function X^(x) reduces to the spherical Bessel 
function, and equation (3.18) reduces to 



QiiVo, k) 
21 + 1 



Jo 



drj 



{a + 3)ei/2|.(^^i/2^(^^ fc))iz(fcr,) + ehx5c{v, k)ji{kr,) 



+ebvi^ + a)Vcir],k)(^-^j^ji-i{krc) - ^[^•?''+i(^''c)) 



a>i) 



(3.20) 



where Tc is the path of a photon, defined as Tc — rjo ~ rj. 

As we are interested in the statistics of the angular fluctuations at a point of an observer Xq = 
x(r7o), we consider the two point angular correlation function (0(770, xq, 7*)6(?7o, xq, 7* )). Expanding 
the fluctuation O at the point Xq by the spherical harmonics. 



8(?7o,xo,7*) = y^a;mYim(r^-y), 



hn 



we can express the ensemble average of the two point angular correlation function as 
(6(770, xo,y)e(77o,xo,y')) ^^Y.^2l + \)CiPi{cose) , 



(3.21) 



(3.22) 



where we defined C; — (af^), and cos 6* = 7* • 7*'. Here the angular power spectrum Cj is given by (e.g., 
Hu & Sugiyama 1995b) 



dkk^M, 



Qiim, k) 



21 + 1 



(3.23) 



where M/ = {P - K) ■ ■ ■ {P - PK)/{P - A')' and P = P + K. 

As we see in equations (3.23) with (3. IS) or ( 3.20| ), the angular power spectrum C/ consists of three 
terms, which are due to the source clustering (5c, the bulk motion of the sources V^, and the gravitational 
potential fluctuations ^. As we will see in the below, the term proportional to 5c in the right-hand side 
in equations ( 3. IS ) and ( 3.20 ) are the dominant contributors to the fluctuations in the CXB (Lahav et al 
1997). Note that the effect of the motion of an observer, i.e. Compton-Getting effect (hereafter C-G 
effect) is not taken into account in this solution. We integrate the above equation numerically. The 
results are described in the next section. 
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IV. FLUCTUATIONS 



In this section we present our results. We here consider three cosmological models, a standard CDM 
model (SCDM), a CDM model with a cosmological constant A (ACDM), and an open CDM model 
(OCDM). As for the initial density power spectrum, we use the Harrison- Zeldovich power spectrum for 
SCDM and ACDM with the COBE normalization (Bunn & White 1997). In the OCDM model, we use 
the spectrum predicted in a simple open inflation model (Yamamoto & Bunn 1996; White & Silk 1996). 
The spectrum is almost the same as the Harrison-Zeldovich one at the scale of the large-scale structure, 
so that our choice of the spectrum does not essentially alter the results. In Appendix C an useful formula 
for the evolution of the CDM perturbation is summarized. 

We have quite a lot of freedom in the parametrization of the X-ray background model. The parameters 
are classified into two kinds. One is the cosmological parameters and the other is the model parameters 
of the X-ray sources. We present how the results depend on the two kinds of parameters. 

A. dependence on the cosmological models 

Before discussing the dependence on the cosmological parameters, we mention the physical contents 
of the fluctuations. Figure 1 shows the angular power spectrum Ci due to the source clustering, bulk 
motion, and gravitational potential fluctuations. It is apparent that the source clustering term is the 
dominant contributor to the large- angle fluctuations. We have checked that the other terms due to the 
gravitational potential and the bulk motion become important only when z^i^ > 0.1. We therefore take 
only the source clustering into account hereafter, focusing on the case Zmin ^ 0.1. 

Now let us discuss the dependence of the fluctuations on the cosmological parameters. Fig. 2 shows 

1/2 

the dependence of the dipole anisotropy C^J^-^ on Hq. As we arc not interested in the spectrum shape 
of the angular fluctuations here, we focus on the amplitude of C;=i. As expected, ACDM and OCDM 
agree with SCDM in the limit of Oq ^ 1- It is shown that the dipole anisotropy is not sensitive to the 
cosmological parameters except for the case of extreme low value of (Lahav et al. 1997). As we see 
in Appendix D, this feature of the dipole anisotropy is attributed to the difference of the shape of the 
power spectrum. Thus Oq dependence can be changed by the normalization scheme of the matter power 
spectrum. The other point to note is that the amplitude of the fluctuations strongly depends on the 
source parameter Zmin- 

B. dependence on the X-ray sources models 

We have three parameters, Zmin, ^max and p, as the X-ray sources model. The evolution of the X-ray 
sources is governed by p and Zmax- From now on, we set the cosmological parameters (ilo = 1, /i = 0.5) 
for SCDM, (f7o 0.3, Oa = 0.7. h = 0.7) for ACDM, and (fio = 0.3, I^a = 0, h = 0.7) for OCDM. 
In Appendix D, we described an analytic calculation for the angular power spectrum C;. The analytic 
calculation is limited only to the case of Oq = 1- However, this is instructive to understand the behaviour 

1 /2 

of the angular power spectrum. Fig. 3 shows the dependence of Cji^^ on the evolution parameter p. 
The dipole anisotropy drastically drops with p. The fluctuation is roughly determined by dividing the 
anisotropic X-ray flux A/^ by the isotropic one li^K As shown in Appendix D, the dependence on p comes 
from the isotropic component li^^ of the CXB. The increase of p, i.e, the increase of the redshift evolution 
of sources leads to the relative increase of /i"-* compared with A7^. Thus the fluctuation decreases with 
p. Fig. 4 shows the dependence of the dipole anisotropy on Zmax- The amplitude of the fluctuation also 
decreases with ^max- The reason is the same as p, but the dependence on Zmax is weaker than that on p. 
These features are common to the higher multipoles (Lahav et al. 1997). 
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Finally, we focus on ^min which was fixed on zero in the paper by Lahav et al. (1997). If the nearby 
bright X-ray sources were removed from the observational data, it would not necessarily be fixed at 
-^min = 0. Fig. 5 shows the dependence of the dipole anisotropy on Zmin- It is apparent that the dipole 
amplitude is sensitive to this parameter especially for Zmin ^0.1. The amplitude of the dipole moment 
rapidly drops with Zjjjin for Znjin ^0.1. In contrast, it becomes almost constant for Znjin ^^ 0.1. This 
strong dependence on Zmin is related with the shape of the matter power spectrum P{k) as discussed in 
Appendix D. The choice of 2;min changes the shape of the angular power spectrum of the fluctuations 
(see Fig. 1). This means that the sources distributed relatively close to the our Galaxy are the dominant 
contributors to the fluctuations in the CXB for the low multipoles. The peak of the matter power 
spectrum is located at A ~ 100ft.~-'-Mpc, which corresponds to z ^ 0.03. The shape of the matter power 
spectrum is the reason why the sources at the low redshift z ^ 0.1 are the dominant contributors to the 
large- angle fluctuations in the CXB. 

V. SUMMARY AND DISCUSSIONS 

In this paper, we have investigated the large-angle fluctuations in the CXB due to the X-ray source 
clustering. We have developed the formalism to describe the CXB fluctuations using Boltzmann equation 
under a simple model of the X-ray sources. Our formalism is a simple extension of that by Lahav et al. 
(1997) to be applicable to an universe with hyperbolic geometry. The dependence of the fluctuations 
on the model parameters has been examined in various cosmological models. The fluctuation does not 
strongly depend on the cosmological parameters. It is quite sensitive to the parameters for the X-ray 
sources (Lahav et al. 1997). The fluctuation is determined by the ratio of anisotropic X-ray flux to the 

isotropic one, i.e., ^/Ui ~ A/y//^"''. A/^ is essentially determined by the nearby sources at low redshift 
for the low multipole moment, while J'-^) is done by the high-z sources. The large redshift evolution of 
the X-ray emissivity (e.g., p = 3, and Zmax is large) makes the flux from the far sources large and, as a 
result, the amplitude of fluctuation becomes small. 

We have also pointed out the importance of the source distribution parameter z^ain- The low multipole 
anisotropics are sensitive to it. This parameter is closely related to the reconstruction of the X-ray map 
by removing nearby bright sources. The dipole moment with Zmin = 0.1 is smaller by order of magnitude 
compared with the case Zmin = 0. If the sources sufficiently close to the our Galaxy 2 ~ 0.1 are removed, 
it may be expected that the effect of the peculiar motion of the observer (Compton-Getting (C-G) effect) 
is dominant in the CXB dipole. 

Let us roughly compare the effects which contribute to the dipole anisotropics except for the source 
clustering effect. The observer's peculiar motion relative to the CMB was measured by using the COBE 
four-year data (Lineweaver et al. 1996), which gave the peculiar velocity Vohs = 368.9 ± 2.5 km/s in the 
direction {I = 264°; b = 48°). Assuming we have a similar motion relative to the CXB, the expected C-G 

dipole is estimated as y^C^ ~ 5.0 x 10~^. 

On the other hand, it is well known that the shot noise fluctuation arises from the discreteness of 
the sources. Since this fluctuation is originated from the Poisson fluctuation of the discrete sources, the 
spectrum is white noise. The shot noise fluctuations in the CXB have been investigated by Lahav et al. 

(1997). Following their result, the amplitude of the fluctuation is estimated ^ Cf^ ~ 1.2 x fm'^, where 

/m is a flux cut-off of bright removed sources in unit of erg s~^ cm~^. If the flux cut-off level becomes 
lower, the amplitude of shot noise decreases. Thus the shot noise fluctuation depends on the flux cut-off 
in observational data and is important when comparing a theoretical model with the observed X-ray map 
(Treyer et al. 1998). The flux cut-off was ~ 3 x 10~^^ erg cm~^ in the HEAO-1, and the amphtude 

of shot noise is estimated as y^C?^^ — 2.8 x 10"'^. Note also that the cut-off level is relevant to the 
parameter Zmin in the sense of bright source removability. 

The dipole owing to the source clustering is shown in Fig. 5 {p = 3, Zmax = 3). In the case of 2;min = 0, 
the dipole anisotropy due to the clustering effect is comparable to the C-G effect. However, if z^i^^ is 
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0(0. 1), the C-G dipole becomes well above the dipole due to the source clustering. If the flux cut-off 
becomes lower by the improvement of observation, the shot noise and the clustering effect could be 
sufSciently smaller than C-G effect to yield a good chance of measuring the C-G dipole in the CXB. This 
investigation also suggests that the careful treatment is required when comparing the observational map 
with the theoretical prediction. When the redshifts of the sources can not be determined, the subtraction 
of nearby X-ray sources from the observed map may contain a delicate problem because the nearby faint 
sources may contribute to the fluctuations sensitively. These problems are left as future problems. 
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APPENDIX A: BOLTZMANN EQUATION 



We write the Boltzmann equation for the distribution function /(t, a;% g, 7*) of photons as 

dt + dx^ dt ^ dqdt^ 97' dt - ^^^''^ ''^'^ ^ ' 

with the source term S. Here we defined the energy in the locally orthonormal frame as q = | | , thus 
{q,q^) forms the 4-momentum of the photon in this frame. Since the Boltzmann equation is written in 
terms of the momentum q measured by an observer in the cosmological rest frame, we must rewrite the 
terms dx^ /dt and dq/dt in terms of and q in order to solve this equation. These relations are given 
from the geodesic equation of a photon. However, the geodesic equation is commonly written in terms 



of the 4-momentum p'^ in the frame (2.1), where is defined by = dx^ /dX with the afhne parameter 



A. The 4-momentum (q, q^) in the rest frame is related to the 4-momcntum p^, as follows: 

9 = (1 + , (A2) 
q' = a(l + $y . (A3) 

The spatial vector 7* is defined as 7* — q^ jq and 7* satisfies "fij^^^^ = 1. 

Equations ([A^) and (^3) give the following relations, up to the first order of and <i>: 



dx^ p^ 1 
'dt 



P (1 + * - $)7' , (A4) 



and 



dq f a* dx^\ r, , ^,dp° 

— = 1 75 + (1 + ^') — . 

dt \dt dxi dt )^ ^ ' dt 

Meanwhile the geodesic equation in the first order of the perturbation is 

dpO (Ida ■ 2 8-9 , ■ Ida 



(A5) 



dt a dt a dx^ ^ ^ a dt ^ ^ 



Inserting this into (A5), we have 



dq / 1 da • 1 9* , , , 



dt \a dt a dx^ 



Thus we can write down the left-hand side of the Boltzmann equation (Al) using equations (A4) and 



( A7). If we employ the conformal time defined by adrj — dt instead of the proper time, it becomes 



^ + (1 ^ $ + *)^7' - + ^ + + = aSit x^ q V) (A8) 



APPENDIX B: MATHEMATICAL FORMULAE IN AN OPEN UNIVERSE 

In this appendix, we summarize the mathematical formulae which are needed in section 3. The results 
of this section are based on the previous work (e.g., Gouda, Sugiyama, & Sasaki 1991; Hu & Sugiyama 
1995b; Wilson 1983). 
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We first consider the scalar harmonics Qk{x) on a hyperbolic universe. It follows (5k|j'' = —k^Qk, 
where \i denotes the covariant derivative in the hyperbolic space with the line element 

j.jdx'dx^ = -^{dx^ + sinh' xdflf^)), (Bl) 

where K = —Hq{1 — Qq — fljy) is the spatial curvature parameter and the line element on an 

unit sphere. Since we assumed that the emitted photons do not scatter, the trajectory of the photon is 
the free streaming. Therefore, the photon geodesies are radial, we may use the radial scalar harmonics 
to describe the free strea ming behaviour. It can be expressed as Qk = Xlj{x)Yim{Q) using the radial 
function X^, in equation ( p.l9| ) and the spherical harmonics Yi^ 



The function Gi for the multipole decomposition in equation ( 3.161 ) is defined 

(0 



G,(x,70 = (-fc)-'Qk|.,....,(x,k)/7--^'(x,y), (B2) 



and 



pii---ii + i _ 2H-1 pi2---ii + i) _ ^ ^Xirh ph---ii + i) 

where the parentheses denote symmetrization about the indices. For the special case of the flat universe, 
it reduces to Gi — (— «)' exp(ix • k)P;(K* • 7*), where Pi{x) is the Legend re po lynomial and = k^/k. 

We derive the evolution equations for the moments Oi in equation ( 3.16 ). For this purpose, we use 
the fact that scalar perturbation quantities are decomposed as ^{ri,x) = $(77, A;)Go(x), etc., and 
JiVc = y^i. Vc iv, fc)Gi(x, 7') for the velocity. The evolution equations for the moments are given from 
equation ( p.l2| ), by the Boltzmann hierarchy: 

eo+e©o--^fcei-(a + 3)$ + e6jf4, (B3) 



3 

Si + (Qi = - 3^Je2 + kOo + (a + 3)k'i' + {a + 3)^bvVc, (B4) 

ei + iei^-^±^k{i~i{i + 2)^)ei+^ + ^^kei.i, for(z>2), (b5) 



where we used 
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APPENDIX C: USEFUL FORMULA FOR EVOLUTION OF CDM PERTURBATIONS 



We here summarize an useful equation for the evolution of the CDM perturbation (Peebles 1980). The 
CDM perturbation Sc{ri,k) obeys 



a • _ 3 ^IqHq 
a 2 a 



(CI) 



where the expansion rate a/a is given by equation (2.2). Note that this equation is independent of the 
scale fc, since we are considering the matter dominant stage. Using equation (2.2), we rewrite the above 
equation as 
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{an„ + a^^K + a^n^) ^ + {^^o + 2a^lK + Sa^^iA) ^ - | = 0. (C2) 
The growing mode solution is written as 



- ^\/% + ^+"A fdaY , (C3) 



which is normahzcd as (5c = i at a <C 1. 



APPENDIX D: ANALYTIC APPROACH 



It is very instructive to evaluate the fluctuations in the CXB in an analytic way. To perform the 
analytic calculation, we conside r th e standard CDM model. In the case of 51o — 1, we have a = {il/fjof' 
with rjQ = 2/ Hq from equation (2.2). Here we use the notations rjuiax/vo = (1 + •Zmin)""'^^^ and ?ymin/f7o = 
(l + Zmax)"'^ ^^- O mitting the terms of the source bulk motion and the gravitational potential fluctuations 
in equation ( 3.20| ), we have 



21 + 1 



->x 



df] eScii], k) ji{k{r]Q - ij)). 



From the definition of e, equation (3.13), we can write 



(Dl) 



QiiVo.k) bxScivo,k) r 



2/ + 1 



E 



dr]i 



Wo 



2{2+a-p) 



3i{k{m - v)), 



(D2) 



with 



2(l+a-p) 



m 



3 + 2a-2p 



(l ^min) ^ ^ (-^ '2^max) 



-3/2-a+p 



(D3) 



where we used the relation Sdrj, k) — (5c(%, k)a{ri). Putting (D2) into ( 3.23 ), we have 



Ci 



2 fb 



dk P{k) 



dx 1 



— 

kriaj 



2(2+Q-p) 



(D4) 



where xi = fc(?7o — rymax), X2 = fc(?7o — '7inin), and we wrote P{k) — (^6c{rio, k)'^). Since k integration is 
effective at krjo ^ 1, we approximate it as 



dx ji (x) 



(D5) 



This approximation works within ^ a few xlO% for ^ p ^ 3 and Zy^in ^ 0.1. Of course this transfor- 
mation is exact when p = 2 + a. 

In the case of Z = 1, a; integration can be done. We get 



Ci^i 



2 fb 



dk P{k) [joixi) -jo{x2) 



(D6) 
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Furthermore we neglect the term of jo (2^2) assuming that k integration is effective at krjo 2> 1. We get 

9 /h \^ f°° 

Ci^io,-(-f] dk P{k) jo{x,f . (D7) 



Assuming z^^in <C 1, we have 770 — 77max — z^in/Ho, then equation (D7) can be written as 

^ nb'^Sj, (--1±^^^-^P—X r dkkT{kf3^{kzr.,n? ■ (D8) 



l-(l+^max)-3/2-"+P 

where k = k/Hg, 5h = 1-94 x 10^^ (Bunn & White 1997), and the matter transfer function T{k) is 
introduced. This equation is meaningful to understand how the dipole depends on the X-ray sources 
parameters. The dependence of the dipole on p and ^max originally comes from the isotropic component 
of the X-ray background radiations oc s(r]Q), whereas z^i^ affects only the fluctuation part in k integration. 

Introdu cing a cut off parameter fc,nax in k integration instead of the transfer function, we approximate 
equation ( p8| ) as 

c ^h\b\ (^^-^^^i-^^^^P^-^ j ' ^'"^^ dfcfcjo(fcz.„i„)^ , (D9) 

where fcmax — fcmax/-ffo ~ O(IO^). Since ^0(2;) oscillates at a; 1, fc integration gives different behaviour 
depending on the value fcmax^^min around unity, as follows: 



fcmax 

2 



dkkjo{kz 



'max/-^ (fcmax^^min ^ 1), 

'" lll(2fcmax^min) +7E]/2z^iu (fcmax^min 



where 7e is the Euler constant. Essentially, fcmax corresponds to the peak of the power spectrum, i.e. the 
scale of galaxy large scale structure k~l^ ~ 10 ~ lO^Mpc. The value of C/=i is therefore very sensitive 
to Zmin around O(10~^). 
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FIGURE CAPTIONS 



Fig. 1 — The angular power spectrum of the fluctuations for various effects. The labels '5\ V, and 
'^f' denote ^/Cl due to the source clustering, the peculiar motion of the sources, and the gravitational 
potential fluctuations, respectively. Here wc set hx = by = I- The cosmological parameters are taken 
as h = 0.5 and = 1 for SCDM model, and h = 0.7 and flo = 0.3 for ACDM model. The cases 
Zmin = 0, 0.02, 0.1 are shown. 

Fig. 2 — fio-dependence of VCi^i/bx for ACDM and OCDM models. The Hubble parameter h = 0.7 is 
taken here. Each panel shows the case ^;min = 0, 0.02, 0.1, respectively. 

Fig. 3 — p-dependence of VCi=i/bx- The cosmological parameters are taken as h = 0.5 and fio = 1 for 
the SCDM model, and as h = 0.7 and Qq = 0.3 for the ACDM and OCDM models. Each panel shows 
the cases ^;min = 0, 0.02, 0.1, respectively. 

Fig. 4 — 2;i„ax-dependence of \fCi^\lbx- The cosmological parameters are the same as Fig. 3. 
Fig. 5 — Zmin-dependence of VCi^i/bx- The cosmological parameters are the same as Fig. 3. 
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